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ABSTRACT 


With the requirement for a spacecraft to maintain an orbital altitude band, a 
simple energy balance algorithm has been developed using a combination of radial 
distance and spacecraft specific energy for fixed-vector thruster control. While each 
trajectory produces a unique band, initial attempts at producing a pre-specified band have 
been unsuccessful. It is theorized that a certain radial bandwidth would correspond to a 
specific set of control parameters, and that by creating maps of the relationship between 
the two for various spacecraft configurations a method of maintaining the pre-specified 
band could be found. This thesis studies variations in spacecraft configurations and finds 
dependence of orbital bandwidth on thrust-to-drag ratio and ballistic coefficient. Also, 
within certain ranges of the control parameters, multiple trajectories produce equivalent 
radial bands. Analysis shows that all single-burn trajectories are characterized by similar 


efficiencies, and are less efficient than a Forced Keplerian Trajectory (FKT). 
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I. INTRODUCTION 


Optimizing orbit maintenance trajectories for low-Earth orbiting (LEO) spacecraft is 
increasingly important as space agencies around the world strive to maintain orbiting 
programs with shrinking budgets. Conceptually, an orbit may be maintained using a 
Forced Keplerian Trajectory (FKT), where thrust is used to cancel drag. This would prevent 
orbital decay, maintaining a true Keplerian orbit. While this strategy is technically 
impractical due to limitations in thrust vectoring and thrust magnitude adjustments, it does 
provide a benchmark by which other strategies may be measured. While optimal control 
theory [Ref. 1], shows that an FKT is not the optimal solution to the orbit maintenance 
problem with respect to fuel, studies have yet to find a trajectory which is more fuel- 
efficient. 

Historically, most orbit reboost strategies have been based on the Hohmann transfer, 
consisting of two thruster burns, where the first burn starts the spacecraft on an elliptical 
path to a higher altitude, and the second bum circularizes the spacecraft's orbit at the desired 
altitude. The Eccentricity Intercept Targeting and Guidance (EITAG) trajectory developed 
by Gottlieb [Ref. 2] notes variations of eccentricity, e, vs semi-major axis, a, and integrates 
backwards in time from the desired orbit. A two bum trajectory is developed which reboosts 
and circularizes the orbit while minimizing errors in achieving the target orbit. 

A third orbit maintenance strategy proceeds from the notion that a spacecraft must be 
maintained within a specific radial band. A space station, for example, must be maintained 
within an altitude band that can be reached by various launch vehicles for servicing and 
manning requirements. A precise, circular orbit is not necessary. With this in mind, a single- 


bum trajectory, boosting the spacecraft as it nears the lower altitude boundary, then allowing 


it to decay over time, would maintain the required band while reducing the number of 
required burns, and perhaps, the fuel required for orbit maintenance. 

Pauls [Ref. 3] and Wilsey [Ref. 4] developed a computer model for such a strategy, 
using various parameter variations to optimize the trajectory. While their results show that a 
single-burn trajectory is less fuel-efficient than the FKT benchmark, further study is 
warranted. In addition, an accurate method of predicting the actual controlled radial band is 
desired. 

By non-dimensionalizing the equations of motion, greater clarity and precision will be 
achieved when varying the control parameters. These results will be used to create maps of 


radius and energy which will allow better control of the spacecraft's radial band. 


II. DEVELOPMENT OF THE ORBIT MAINTENANCE PROCEDURE 


A. NON-OPTIMALITY OF A FORCED KEPLERIAN TRAJECTORY 

While optimal control theory techniques show that a Forced Keplerian Trajectory 
(FKT) is not optimal with respect to fuel [Ref. 1], significant insight into the problem can be 
gained by looking into an extremely simplified version of the orbit maintenance problem 


[Ref. 5]. If an FKT is optimal, the theory should lead to an optimal control history where 


thrust is equal to drag. 





Figure 1 Orbital band 


The problem may be simplified by stating that the only requirement of the control system is 
to maintain the orbit semi-major axis, a, and eccentricity, e, within such limits that the orbit 
remains within the specified band. These two orbital elements define the shape of the orbit 


within the orbital plane. This shape could theoretically vary from a circular orbit (e=0) 


where a varies from Rin to Rmax (see Figure 1), to an elliptical orbit with perigee equal to 
Rin and apogee equal to Rmax , with e variation dependent upon the size of the band. 

Further simplification is possible by looking only at the portion of the problem relating 
to semi-major axis, a. This leaves a one degree-of-freedom system for optimization. To get 
an equation of motion for semi-major axis, the work-energy theorem is used 


AE= [F-ds (1) 


This can be used to relate the specific mechanical energy of the orbit to the forces involved, 


namely thrust and drag 
de _(T-D ds 
dt m )dt 


=) (2) 





Recalling the definition of the specific mechanical energy, © = - / 2a, 








de uw da 
a 3 
dt 2g? dt (3) 
This leads to an equation for a, 
a=22%(T_p) (4) 
pm 
The other equation of motion for the system comes from the rocket equation 
vee ee ae (5) 
LB. ae 


where the specific impulse, I.) and the acceleration due to gravity are multiplied to define the 
exhaust velocity of the engine, Ve. 

The final major simplification comes from focusing on the band to be maintained. 
When attempting to maintain a small enough band, the change in a is small, so the squared 
term is essentially constant. Additionally, such a band would require a small change in 
velocity, Av, to maintain it. Thus, v is also assumed to be constant. These assumptions allow 


Equation 4 to be rewritten as 


ae“ (T-p) (6) 
m 


where k is a constant. 

The optimal control problem is now formulated as finding the control history that best 
transfers the spacecraft from some initial condition a, to some final condition ay. The 
optimality criterion is minimizing fuel burned, which can be stated as maximizing the final 
mass. The problem constraint is the available thrust, which is defined as 


O<T<T., (7) 


Thrust is the single element of control in the simplified problem. From optimal control 


theory, the Hamiltonian is 


2 T- + (8) 


where A, and A. are the Lagrange multipliers, or costates of the system, corresponding to 


the semi-major axis and mass respectively. 
Pontryagin's Maximum Principle (PMP) states that at every point on the trajectory of 
the system, the optimum control history maximizes the Hamiltonian. Stated mathematically, 


u’ =argmax, _,H (9) 


By examining Equation 8, it is clear that the control variable, T, that maximizes the H- 
function depends upon the sign to the quantity in brackets, which is called the switching 


function, s 


ACKOK 
§ 2 ee on (10) 
mY; 


Noting the constraint upon thrust, application of the PMP shows that 


0 
re = (uy) 
0 s<@ 


Called a bang-bang controller, this optimal control history contains no reference to drag; 
therefore the optimal control for such an orbital transfer does not define thrust equal to 


drag. The Forced Keplerian Trajectory is thus shown not be the optimal control method. 


B. ORBIT CONTROL STRATEGY 

It is desired to maintain a spacecraft within a specified radial band, using a single 
thruster firing to reboost the spacecraft when it nears the minimum radius. Wilsey [Ref. 4], 
shows that such a scheme is possible using radius and spacecraft specific energy as control 


parameters. Figure 1 shows such a radial band, defining radii and energies. 


To facilitate this orbital control scheme, it is necessary to initiate thruster firing prior to 
the spacecraft decaying to the minimum altitude, and shut down thrusters prior to gaining 
the maximum altitude. Thruster firing is controlled by orbit radius, with 

ér=R,-1 (12) 


th 


where ryp is the radius of thruster initiation. Spacecraft specific energy is then monitored 
throughout the burn, with thruster cut-off occurring when 

E=En (13) 
with E,,, being specific energy at thruster cut-off. By specifying values for dr and Exyy, 
minimum and maximum radii will be achieved, denoted by rmy and rmx, respectively. This is 


pictured in Figure 2. 





Figure 2 Thrust Control Parameters 


The choices for dr and F,;, will produce a unique maximum band, but what values 
correspond to the desired band? Does a specific, desired band produce a specific trajectory? 
How do the bands and associated trajectories vary with changes in spacecraft parameters 
such as thrust-to-drag ratio. specific impulse. and ballistic coefficient? Ross, Pauls, and 


Wilsey [Ref. 6] propose a simple "energy-balance" algorithm where Or and Fy; are varied 


until the desired radial band is achieved [Appendix. A]. By following this procedure and 
mapping Or and E,p as they vary with the actual band attained for different spacecraft 
configurations, a relationship between the parameters is sought which would remove the 


need for iteration and allow orbit trajectory prediction. 


C. DEVELOPMENT OF THE EQUATIONS OF MOTION 

The initial investigation of a single burn, low-earth orbit trajectory will be simplified by 
assuming an orbit which is initially circular. The spacecraft is a non-lifting (blunt) body, 
therefore the only aerodynamic force which affects it is drag. The external forces thus 
perturbing the orbit are aerodynamic drag, gravity, and thrust. Looking at the orbital 
dynamics as a two-body problem, the coordinate system may be defined as shown in 


Figure 3. 


transverse 





The equations of motion for this system are 


F 
ay ie (14) 


where a, and a,, are the radial and transverse components of the inertial acceleration, IF, 
and ZF, are the sums of the respective forces, and m is the spacecraft mass. As seen in 
Figure 3, drag and thrust can each be broken into components 

D,=-Dsin(y} D,=-Dcos(y} (16a&b) 


T,=Tsinja) T,=Tcos(a} (17a&b) 


where y denotes the flight path angle (between the velocity vector and the transverse axis), 
and a represents the thruster elevation angle (between the thrust vector and the transverse 


axis). Rewriting the equations of motion in polar coordinates using equations 16 and 17 





leads to 
ea p D T 
#-Ap=-— 14 (18) 
rom m 
4 : De iT 
@r+ 267=-_*+ 2 (19) 
m m 


where w is the earth's gravitational parameter, and 6 is the angular position of the spacecraft. 


D. NONDIMENSIONALIZATION OF THE EQUATIONS OF MOTION 

When dealing with an orbital system, the use of standard dimensional units leads to 
large variations in numbers that cause computational errors. To counter this and to better 
study the effects of variation of certain parameters, a system of canonical units is developed 


to non-dimensionalize the equations of motion. 


1. Definitions 
To nondimensionalize the equations, base units for length, mass, and time are 
defined. The base length is the radius of the target orbit, which is defined as the midpoint of 
the desired radial band, and is denoted by rp. The base mass is the initial spacecraft mass, 
enoted by mp. The base time is defined as the period of a circular orbit at a radius of 
Tp divided by 2x, and is denoted by ty, 
r 


t= re (20) 


These base units are then used to define the nondimensionalized variables for radius, mass, 


and time 


y (21) 
m= 
me (22) 
Tat 
‘p (23) 


The variable 6 is dimensionless by definition. 
In addition to the three "basic" units, a base unit for density is also defined, and denoted 


by pp. Using a local exponential density model with atmospheric scale height B 


eho 


P=P, (24) 
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Pp is defined as equal to p,. Nondimensionalizing B with rp . and setting rep equal to rp 


produces a nondimensionalized density 
pwe ty (25a) 
where 
B=Br, (25b) 


This is used in the determination of drag, as drag is given by 
(26) 


1 = 
D=—pAC,v 
where A is the surface area, Cg the drag coefficient, and v the spacecraft velocity. A table of 
quantities, their nondimensionalized cousins and relating factors is given in Appendix B. 


2. Equation Nondimensionalization 
The above definitions are used to develop the nondimensionalized variable 


derivatives 
r wt 
(27) 





b--F(S eS Bus ape) oN: A (28) 
ala) ala, ade, 1 & 


where primes denote a nondimensionalized parameter differentiated with respect to time 


These derivatives are used to yield the nondimensionalized equations of motion 





T 
- ot (29) 
rom sm 
7 7 D, T 
O'T+ 207 =-—1 += 
m m 


(30) 


where the nondimensionalized values of both drag and thrust components are related to their 
dimensionalized counterparts by the same factor 


= t° 
F=F2U F=DT 
mf 


bob 





(31) 
3. Parameter Development 


To study the effects of varying spacecraft design parameters on the orbit trajectory 
and their related bands, nondimensionalized parameters are developed which represent 
ballistic coefficient, engine specific impulse, and thrust-to-drag ratio. These parameters can 


be thought of as quantifying atmospheric effects, spacecraft shape, and engine sizing. 
a. Nondimensionalized Ballistic Coefficient, B 
Recalling equation 26, the nondimensionaliztion of drag requires the 


nondimensionalized density defined above, p , anda nondimensionalized velocity, Vv 
Using the velocity equation 


2.2 242? 
ver +r 80 


and the relationships in equations 27 and 28 gives 


(32) 
vert (33) 
t 
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which is substituted, along with p into equation 26, Nondimensionalizing drag using 
equation 31 gives 


(34) 





By defining the ballistic coefficient as the base mass divided by the product of the surface 


area and the drag coefficient, 


Bev (35a) 


the nondimensionalized ballistic coefficient is then defined as 











pe (35b) 

Pt, 
which leads to the final nondimensionalized equation for drag 

pet P* (36) 

“ B 
b. Mass/Thrust Parameter, p 
Fuel mass flow and thrust are related by the rocket equation, 
ieee (37) 
Is 
spo 


Following the example of equation 27 when nondimensionalizing the derivative of mass and 


using the thrust relationship of equation 31, the nondimensionalized mass flow is written 


13 








i= -T—— (38) 
re 
By defining 
T 
p-— (39) 
1 8ot, 


variations in spacecraft specific impulse, I, , can be incorporated into a final 
nondimensionalized equation for both mass flow and thrust 


in’ =— pT (40) 


c. Thrust-to-Drag Ratio, T, 


The thrust-to-drag ratio is defined as the thrust divided by the drag at the target 
orbit. As both thrust and drag are forces and thus nondimensionalized by the same factors, 


this ratio may be written 


(41) 


4 
u 
o| ji 


IH. COMPUTER MODEL DEVELOPMENT 


A. ORIGINAL DEVELOPMENT 

Originally developed by Pauls [Ref. 3] and refined by Wilsey [Ref. 4], a computer 
program, written in FORTRAN, was developed to simulate the orbital motion of a 
spacecraft. The main routine controlled input, output, and subroutine calls. Subroutines 
computed drag and orbital parameters, controlled thrust, and updated the equations of 
motion. A fourth-order Runge-Kutta numerical integration routine was used to integrate the 
equations of motion. The nondimensionalized equations were simplified by defining new 
state variables, which transformed the two second order differential equations into a set of 
four first order differentials. 
B. STATE VARIABLE DEFINITIONS 


A five element state vector, x , is defined as 


X,=r (42) 
x,ote (43) 
x,=6 (44) 
x= 0-2 (45) 
x,=m (46) 
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t 








These states can be nondimensionalized using the factors introduced previously and 
substituted into the equations of motion. Combining equations 29 and 30 with the 
nondimensionalized state variables produces 
dx, 4 Ds. T 
Sax 7-1 -Psinys Tsing (47) 
dt x om m 
dk 2% ,X, D 
—=-_ - _D_ cosy + ——— cosa (48) 
dt x x,m X 
As the flight path angle, y, is a function of the velocity vector, the geometric functions can be 
represented by 
sin y=— cos y=—L— (49a&b) 
v ¥ 


Substituting these relationships into equations 47 and 48 leads to the final equations of 
(50) 


motion as they are programmed into the simulation: 
dk, 
= x, 
(51) 














dK 250 4 
——=xX xX,-— — 
df x° mv 
dk, (52) 
—e = X 
dt ° 
dk 2x,%, Dx T 
bacon Xa%e - ne + ate cosa (53) 
x mv xX,m 


| (84) 


C. PROGRAM EXECUTION 

Program execution begins with the input of variables from the user as well as an input 
file (see Appendix A). Following the definition of the initial conditions, the DRAG 
subroutine is called. The first call establishes the drag of the FKT benchmark. Execution 
then enters the main loop, where drag is called again, computing the instantaneous velocity 
as well as the density. With the velocity updated, the specific energy is computed and the 
CONTROL subroutine is called. This routine controls thruster activation and cut-off. With 
a thrust value thus determined, subroutine EQN updates the equations of motion. RK4 is 
the fourth order Runge-Kutta routine that integrates the state variables. Flight-path angle, 
fuel-burned totals and achieved altitudes are updated once every radian, while the routine 
ORBELMTS updates the classical orbital elements. The equations of motion are updated 


5000 times per orbit, while all pertinent variables are sampled ten times per orbit. 


D. VALIDATION 

Since the program logic remains intact, the validation process is simplified. The initial 
baseline test consists of running the program with the spacecraft in an initially circular orbit 
with gravity the only external force. Radius, velocity, specific energy and angular 
momentum all remain constant, validating the basic equations without drag or thrust. 

The next step in the validation process compares the trajectory produced by the 
modified version with those produced by the earlier program for a consistent set of inputs. 
Since an in-depth validation is done by Wilsey [Ref. 4], duplicating the results of the earlier 
version validates the modified program. The main difference between the two versions lay in 


the nondimensionalizing process, and the calculating of instantaneous spacecraft mass. To 
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achieve universality, all input variables and parameters are already nondimensionalized. 
Using input values which correspond to those used in the earlier study, the modified version 


produces equivalent results. 


E. DEVELOPMENT OF "SMART" VERSION 


The aim of this study is to determine values of dr and E,, which will maintain a 
specified band, therefore a version of the program that accepts input values of dr and E,, to 


simulate the trajectory is developed. This program is used as a test to validate the results 
obtained using the parameter mapping procedure. A listing of the program is found in 


Appendix C, while example input and output files are listed in Appendix D. 


IV. PARAMETER MAPPING 


Toremove the iteration on 6r and E,, which is required in the algorithm, a way of 


determining these control parameters for a given band prior to program execution is 
required. In addition, the variations in the control parameters for different spacecraft 
configurations must be known. 

Toachieve these goals and possibly find a mathematical relationship between the 
parameters, a mapping technique is used. Running the simulation for a specific spacecraft 
configuration produces a set of trajectories. Multiple simulation runs for various 
configurations produces a family of trajectory curves. Combining all of the trajectories into 


a catalog where the achieved band is graphed against 6r and E,, would bypass the 


algorithm's iteration process; for a given spacecraft configuration and a desired orbital band, 


values for dr and E,,, could be pulled from the maps and input into the control program. 


A. PARAMETER VARIATION 

For this work, a spacecraft configuration is defined as a set of the three 
nondimensionalized parameters previously introduced. Five values for each parameter are 
combined to produce a database of 125 different configurations. The parameter values and 
their dimensionalized counterparts are listed in Table |. These values are chosen to 
represent a range of actual values for thrust-to-drag ratio, ballistic coefficient and specific 
impulse (Ip). Increasing values of T/D represent moving up to more powerful motors. 
Increases in specific impulse, the amount of thrust per unit weight of propellant, are 


reflected in decreasing values of p. The I, values reflect those of different propellent 


systems, from cold gas propellants (1,, = 100), to ionic propulsion systems (Isp = 3000). 
Increasing values for B , the nondimensionalized ballistic coefficient, represent an increase in 


Table 1 PARAMETER VALUES FOR MAPPING PROCEDURE 
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spacecraft mass, decrease in surface area, or both, according to equations 35a and b. A 
representative large orbiting platform with the following characteristics 
mass = 20,000 kg 
surface area = 60 m2 
drag coefficient = 2.2 
would have a ballistic coefficient of 150 kg/m2. 
B. PROGRAM MODIFICATION 


Two modified versions of the computer simulation are used to accomplish the mapping. 


One tracks the nondimensionalized orbital decay, df, while the other focuses on the thruster 


cut-off energy, E,,. 


1. Automatic variation of ordered band 


Both control parameters, 6f and ae can be defined by an “ordered” band. The 


amount of orbital decay defines the lower bound, while the thruster cut-off energy defines 
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the upper bound. By varying this ordered band, a range of values for 6r and E. can be 


simulated and mapped. A program loop added outside of al! subroutine calls and data 
output commands automated repetition of the simulation, allowing the ordered band to vary 
from SO km to 0 km (which would simulate the FKT), in 2.5 km increments. For each 
ordered band, the program simulated 50 orbits, or approximately three days of spacecraft 
flight. 

2. "Bubble sort" of minimum and maximum radii 

The actual band of a trajectory is defined as the radial distance between the 

minimum and maximum radii that is achieved during the simulation. Two 


nondimensionalized variables,t_ | and Tare initially defined as equal to the initial orbit 


radius. As the program simulates the trajectory and the instantaneous orbit radius is 
computed, it is compared to the two variables. Whenever the instantaneous radius is less 


than the radius heid inf __, the variable takes on this new, lesser radius. Conversely, F_. is 


mn’ 
updated to the instantaneous radius whenever the orbit goes above the current maximum. At 
the end of a trajectory simulation, fF, and Twill hold the minimum and maximum radii 
achieved. The variables are reset for each trajectory, producing an actual band corresponding 


to each variation in 6f and E., 


3. Actual band vs. ordered band 


The purpose of the "ordered" band is tutorial; it provides a way to visualize and vary 


the control parameters 6r and Ee Once the simulation is run, the ordered band has no real 


significance. The actual radial band that is achieved for the given set of control parameters is 
all that is of practical interest. An orbit maintenance problem as developed here is defined by 


the band to be maintained, not the individual obit radii. Any specified minimum and 
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maximum orbit radii form a radial band. By looking at the achieved band rather than the 


minimum and maximum radii separately, the dimension of the problem is reduced. 


C. MAPPING ACTUAL BAND VS. dr AND E é 


Output from each of the parameter-mapping versions consists of the values of the 
respective control parameter and the corresponding actual band produced. An example plot 
for thruster cut-off energy is shown in Figure 4. This depiction of the trajectory is then 
combined with similar trajectory curves corresponding to different spacecraft 


configurations, producing a catalog of trajectories. With a specific spacecraft and a required 


band to be maintained, values of the control parameters of and E. can be pulled from the 


plots and used to maintain the spacecraft’s orbit, removing the need for iteration in the 


proposed al gorithm. 


Achieved Band 


Figure 4 Example map of achieved band vs. thruster 
cut-off energy 
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Vv. ANALYSIS AND RESULTS 


Analysis of the parameter maps consists of investigating the effects of spacecraft 
parameter variation on achieved band. While 125 different spacecraft configurations are 
catalogued, investigation is focused on four configurations which best exhibit the noted 
phenomena. The following configurations are studied: 


Config 533 Config 234 Config 342 Config 453 
T/D = 500 T/D=75 T/D = 100 T/D = 250 
B= 150kg/m2) B= 150kg/m2 B=300kg/m2 B=S5S00 kg/m? 
Isp = 500 sec Ip = 1000 sec 1,5 = 250 sec Isp = 500 sec 





The prefix D- or E- on the maps indicate the source file of the data; D represents actual 


band plotted against orbital decay ( 6r), while E represents actual band plotted against 


thruster cut-off energy (E,,) The first digit refers to T/D ratio, the second to B, and the third 


to I,y. On each plot, only one parameter is varied. 
Accuracy of the mapping technique is also tested by using the control parameters in the 
“smart” version of the simulation and noting the size of the actual band achieved. All figures 


discussed are located in Appendix E. 
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A. VARIATION OF SPECIFIC IMPULSE, I,p 
Figures | through 4 show the effects of varying a spacecraft engine's specific impulse. 
In all cases, the effect of varying I,p is extremely small. Figure 1 shows no variation between 


trajectories. It is also noted that the curves for actual band vs. 6r exactly match the shape of 


the actual band vs. E,, . This is due to the way the control parameters are varied together, 


using the concept of the "ordered" band introduced in the previous chapter. 

Figures 2, 3, and 4 all show families of trajectory curves corresponding to spacecraft 
whose configurations differ only in the value of their engine's specific impulse; thrust-to- 
drag ratio and ballistic coefficient remain constant within each figure. While individual 
trajectories differ, the general shape in each case remains the same. Also, as the actual band 
narrows, the trajectories become practically identical. This convergence represents identical 
control parameters producing equivalent bands for different spacecraft configurations. 


Recalling the definition of the mass/thrust parameter p (equation 39) and Table | it is 
clear that despite large changes in specific impulse, the parameter p remains a smal! number 


that varies little. Figures | through 4 confirm that variations in specific impulse have little 


effect on the orbital band maintained using the single-burn algorithm. 


B. VARIATION OF BALLISTIC COEFFICIENT 

Figures 5 through 8 show the effects of varying the spacecraft's ballistic coefficient. 
Again, the values of thrust-to-drag ratio and (in these cases) specific impulse remain 
constant throughout each figure. Comparing the trajectories in any of the individual maps 
shows considerably more effect on actual band by variation of ballistic coefficient than was 
produced by varying specific impulse. 

Generally, the size of the actual band for a given control parameter decreases with 
increasing values of the nondimensionalized ballistic coefficient. Figure 5 shows this trend 
with only small deviations. Figures 6 and 7 show larger variations between individual 
trajectories, but the shrinking of the actual band does continue. Figure 8 shows less 


variation, which points towards a relationship between the ballistic coefficient and the thrust- 
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to-drag ratio. Varying ballistic coefficient has less of an influence on configurations with 
larger values of T/D. 

Recalling equations 35 and 36, increasing spacecraft mass or reducing spacecraft 
surface area leads to an increase in ballistic coefficient, which in turn decreases drag. As the 
thrusters modeled in this simulation operate at a constant thrust, those configurations with 
large thruster values feel little effect from the decrease in drag, while those with smaller 
thrusters, and thus smaller T/D ratios, more easily “sense” the decreased drag. With thrust 


and drag smaller, the algorithm is better able to control the actual band. 


C. VARIATION OF THRUST-TO-DRAG RATIO 

Figures 9 through !2 show the effects of varying the spacecraft's thrust-to-drag ratio. 
Each plot represents a family of spacecraft with consistent values of ballistic coefficient and 
specific impulse, but different T/D ratios. Comparing individual trajectory curves within a 
map shows that increasing thruster size generally increases the actual orbital band achieved. 

As with ballistic coefficient variation, the effects of varying T/D are more profound than 
when changing specific impulse. Once again, the effect of large thrusters is most noticeable, 
as they seem to cancel out irregularities present in configurations possessing smaller 
thrusters. In each figure, the trajectory corresponding to that configuration with the largest 
thrust-to-drag ratio ( T/D = 500) is a smooth curve, and in some cases almost a straight line. 
Larger values for T/D ratio inhibit the single-burn algorithm; for a given value of either 
control parameter, a larger thruster produces a larger band. Smaller thrusters, which may be 


thought of as being more responsive to this control scheme, produce a tighter orbital band. 


D. MULTIPLE TRAJECTORIES FOR A SPECIFIED BAND 
An assumption made during the formulation of the single-burn algorithm and the 


mapping procedure covered here deals with the relationship between the control parameters 
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and the band they produce. Given a spacecraft and a set, periodic thruster firing schedule 


(defined values of 6f and E. ), a certain radial band will be maintained. The parameter 


mapping technique essentially reverses this process, starting with the required radial band 
and looking at the required thrusting scheme. 
It seems reasonable to assume that since a unique thruster firing schedule produces a 


specific radial band, the converse will also hold; a certain radial band will yield unique 


values of or and E.,- As can be seen in Figures 2 through 12, this is not the case. On each 


map, there is a region where a specific band can be produced by multiple values of each 
control variable. Such a choice between control parameter values could lead to new 
optimizing opportunities: which trajectory is “best"? To investigate this development, the 
"smart" version of the simulation is used to study different thruster firing schedules for a 


specific spacecraft configuration. 


E. "SMART" PROGRAM VALIDATION 

Figure 13 shows the two plots relating actual band to the control parameters for 
Configuration 234 introduced previously. For nondimensionalized radial bands between 
approximately .0076 and .0095 times the base orbit radius, there are three possible values 
for both control parameters. To validate the smart program and investigate the control 
parameter relationships, the following trajectories are chosen which produce nearly equal 


radial bands: 


Trajectory Actual band 


Traj.1 0.008610134 
Traj.2 0.008608818 
Traj3 0.008789484 
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The values of SF and Ex corresponding to each trajectory are used by the smart program to 


generate the specific simulation for 50 orbits. 


1. Trajectory comparison 
The orbit radius of each individual trajectory is shown in Figures 14 through 16. A 
combined plot is shown in Figure 17. While each trajectory produces a radial band of 
comparable width, the bands and the thruster firing patterns are quite different. Figure 14 
shows Trajectory I's decay to initial thruster firing taking just over ten orbits, followed by a 
long burn. After thruster cut-off the orbit's eccentricity is quite significant, as seen in the 


rhythmic variation in nondimensionalized radius (Rbar) as the orbit decays again tof, . 


Figure 15 shows a similar pattern for Trajectory 2, with a smaller orbital decay and a shorter 
burn-time. 

Figure 16 shows a much different pattern. Decay to the initial thruster firing occurs 
rather quickly, with the orbit proceeding directly toa highly eccentric state. The eccentricity 
of the orbit causes the spacecraft to hit the specified thrusting radius as well as the thruster 
cut-off energy more frequently. The orbit actually becomes an ellipse whose apogee and 
perigee are bounded by the minimum and maximum radii of the orbital band. Noting the 
minimum and maximum radius values for Trajectory 3 plotted in Figure 16 and the 


definition of eccentricity 


e= mx mn (55) 


the eccentricity is found to be approximately 0.0043. Figure 18 compares the eccentricities 
of the three trajectories, and that of Trajectory 3 confirms the computed result, as it tends 


towards a steady-state value of approximately 0.0045. 
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2. Efficiency comparison 

Given these differing burn patterns, which one is best? A performance index must 
be defined to be optimized; this study looks at maximizing fuel efficiency. Figure 19 shows 
the mass of fuel burned for each trajectory, as well as that for a Forced Keplerian 
Trajectory. 

Itis clear that the conclusions made by Wilsey [Ref. 4] are confirmed; none of the 
single-burn trajectories are more fuel-efficient than an FKT. What is of interest, however, is 
that the average slope of the fuel-burned curve appears independent of thruster burn pattern. 
Regardless of the frequency and length of thruster firing, all three single-burn algorithm 


produce the same average fuel burn rate. 


3. Variation Of Orbital Simulation Length 
In comparing these trajectories, it is helpful to look at different orbital parameters, 


including spacecraft specific energy. The nondimensionalized specific energy is defined as 


-_ vd 
2. °F (69) 


and is shown for each orbit in Figure 20. The near vertical increases in specific energy 
correspond to thruster firings; the greater the rise in energy, the longer the burn. Trajectory 
3 quickly reaches a state where the specific energy is tightly controlled by a high-frequency, 
short-burn schedule. With only two burns showing for the other trajectories, however, 
sound conclusions cannot be drawn; more data is required. An additional series of 
simulations is run for a period of 100 orbits, with orbit radius, specific energy, and fuel 


mass burned shown in Figures 21 through 24, respectively. 


a. Transition to a Tight Energy State 

From Figure 21, it is seen that Trajectory 2 undergoes a significant change around 
its 65th orbit; the radial band maintained expands significantly. Comparing the three 
trajectories in Figure 22, it is seen that Trajectory 2 transitions to a tightly governed energy 
state very similar to that reached earlier by Trajectory 3 (although the energy level of 
Trajectory 2 is higher). This transition appears due to the eccentricity of the orbit increasing 
to a point where the orbit grazes the prescribed band; the perigee and the apogee of the orbit 
are consistently in contact with the bounds defined by the control parameters, as shown in 
Figure 23. Despite the change in the thrusting pattern and the increase in radial band size, 
Figure 24 shows that the slope of the fuel-burned curve remains consistent, above that of the 


FKT. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The purpose of this thesis was to create a catalog of control parameter maps for various 
spacecraft configurations, using a proposed single-burn trajectory algorithm for low-earth 
orbit maintenance. Such maps are used to remove the iteration required in the algorithm 
(Appendix A), assuring maintenance of the desired band and removing the need for human 
intervention in the control process. Additionally, study of the relationships between control 
and spacecraft configuration parameters led to new insights into the orbit maintenance 
problem. Finally, comparison of all single-burn trajectories to the Forced Keplerian 
Trajectory (FKT) benchmark for optimization of fuel-efficiency could not corroborate the 
optimal control theory finding that the FKT is not the optimal trajectory. 

The FORTRAN simulation developed by Pauls [Ref.3] and refined by Wilsey [Ref. 4] 
was further modified by nondimensionalizing the entire simulation, then developing two 
versions to accomplish the parameter mapping. A "smart" version was also developed which 
would accept as input the control parameters mapped out in the thesis, thus testing the 
validity of the plots. 

Three nondimensionalized spacecraft configuration parameters related to thrust-to-drag 
ratio, ballistic coefficient, and specific impulse were defined and varied through a range 
practical values. Their individual effects upon the spacecraft's trajectory were studied. 
Variations in specific impulse produced little change in the trajectories. Variations in 
ballistic coefficient and thrust-to-drag ratio produced more pronounced changes between 
trajectories, with a coupling noted between the two. While increasing the value of the 
ballistic coefficient generally decreased the size of the maintained band, the influence of the 


ballistic coefficient was significantly reduced for high thrust values. High values of T/D 
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produced larger radial bands, and also seemed to cancel out any orbital irregularities present 
in spacecraft with smaller, less powerful thrusters. 

An interesting discovery was that a specific radial band could be produced by different 
combinations of control parameters; the relationship between an orbital band and a thruster 
firing schedule was not one-to-one. This allowed for further optimization studies, as a 
certain single-burn scheme might actually be superior to an FKT. 

While the FKT was shown to be consistently better than all single-burn trajectories in 
this study, the differences in thruster firing patterns proved intriguing. Different trajectories 
producing comparable radial bands had drastically different thruster firing times and 
frequencies. High orbit eccentricity, where the orbit's perigee and apogee were bounded by 
the radial band , sent the trajectory into a periodic, high-frequency, low burn-time thrusting 
scheme, where specific energy varied little. Fuel-burn curves for each single-burn trajectory 
maintained a consistent slope, even after the orbit transitioned to this tightly governed 
specific energy state. 

This thesis opens the door for further study as well, since some questions were 
unanswered, and new questions have arisen. A comparison of the single-burn algorithm toa 
dual-burn algorithm such as EITAG [Ref. 2] is warranted. While not superior to an FKT, a 
single-burn orbit maintenance algorithm may still be better than the dual-burn solution. 
Additionally, what effects do other orbital parameters such as flight-path angle and 
spacecraft velocity have upon the trajectory? Is the transition to the tight energy state 
globally valid? And ona grander scale, what is the optimal orbit maintenance strategy? What 
form of bang-bang control will minimize fuel burned and provide the most efficient means 
of maintaining a continued presence on orbit? With man's continued presence in space 
seemingly contingent upon maximizing performance and efficiency while minimizing cost, 


such study is certainly warranted. 
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APPENDIX A 


Single-burn LEO Maintenance Algorithm 


_ guess input IC's 
drand Eth and BAND 
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APPENDIX B 


Nondimensionalization Table 
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- PROGRAM LISTING 
PROGRAM ORBMAINT 


APPENDIX C 


ORBMAINT is a self-contained, nondimensionalized orbit maintenance program 
that will control the orbit of one body orbiting another. The program consists of a main 
routine controlling input, output, and subroutine calls, and five subroutines which 
compute atmospheric drag, the equations of motion, the standard orbital elements, and 
the number and duration of thruster firings. A 4th order Runge-Kutta integrator is used 
to solve the equations of motion. The program also performs comparative calculations 
for a Forced Keplerin Trajectory (FKT) 


VARIABLE DEFINITIONS 


ARRAYS: 
XBAR(1) 
XBAR(2) 
XBARG) 
XBAR(4) 
XBAR(5) 
XBDOT(1) 
XBDOT(2) 
XBDOT(@) 
XBDOT(4) 
XBDOT(S5) 


PBAR 
USER INPUT: 


ALPHA 
ALPHAR 
DLTRBR 
ETHBR 
COUNTERS: 
TBAR 
INDEX 
KOUNT 


nondimensionalized orbit radius 
nondimensionalized orbital radial velocity 
theta 

nondimensionalized orbital angular velocity 
nondimensionalized spacecraft mass 
derivative of XBAR(1) 

derivative of XBAR(2) 

derivative of XBAR(3) 

derivative of XBAR(4) 

derivative of XBAR(5) 


initial nondimensionalized orbit radius 
initial nondimensionalized spacecradft mas 
number of orbits to simulate 

time increment step size 

print interval step size 

atmospheric scale height 
nondimensionalized thrust-to-drag ratio 
nondimensionalized ballistic coefficient 
nondimensionalized mass/thrust parameter 


thrust angle 
nondimensionalized orbital decay: ROBAR-RTHBR 
nondimensionalized spacecraft energy 


current time (initial value = 0) 


(initial value = 0) 
(initial value = 0) 
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DOBAR 
EOBAR 
VOBAR 
RMNBR 
RMXBR 
PROGRAM VARIABLES: 
ABAR 
APOGEE 
DBAR 
EBAR 
ECC 
GAMMAD 
GAMMAR 
HBAR 
FUELKEP 
TOTKEP 
MBRFUEL 
PERIGEE 
PRDBAR 
RHOBAR 
RTHBR 
TBARF 
THBAR 
THBRCAP 
THBRKEP 
THGMA 
VBAR 
CONSTANTS: 
PI 
J 


AAA AN AMAANAIAANAANAANANAAANNANNANAANNA 


variable declarations: 


INITIAL AND BOUNDARY CONDITIONS: 


initial nondimensionalized atmospheric drag 
initial nondimensionalized spacecraft energy 
initial nondimensionalizedspacecraft velocity 
nondimensionalized minimum orbit radius 
nondimensionalized maximum orbit radius 


nondimensionalized semi-major axis 
nondimensionalized orbit apogee 
nondimensionalized atmospheric drag force 
nondimensionalized spacecraft energy 
orbital eccentricity 


flight path angle 

nondimensionalized angular momentum 
keplerian fuel bummed per program iteration 
total keplerian fuel burned 
nondimensionalized mass of fuel burned 
nondimensionalized orbit perigee 
nondimensionalized orbit period 
nondimensionalized atmospheric density 
nondimensionalized radius of thruster firing 
simulation stop time 

nondimensionalized spacecraft thrust 
nondimensionalized spacecraft thrust capability 
nondimensionalized keplerian thrust 

flight path angle at thruster initiation 
nondimensionalized spacecraft velocity 


pl 
state variable index 


IMPLICIT REAL*8(A-H,M-Z) 
DIMENSION XBAR(5),XBDOT(5) 


Cc 

C constant definitions: 
PI=DATAN(1!.0D+00)*4.0D+00 
J=5 

C 

C _ initialize keplerian fuel totals 


FUELKEP=0.0D+00 

TOTKEP=0.0D+00 
cceccececccccccccccccceccec 
Cc 


OPEN(10, FILE="input’, STATUS=’"OLD') 
OPEN(20, FILE='orbpar’, STATUS='NEW') 
OPEN(30, FILE="orbel’, STATUS='NEW') 
OPEN(40, FILE="orbetc', STATUS='NEW') 
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C read input file 
READ( 10,1) ROBAR,MOBAR,FORBIT,STEP,PRNT,BETABAR, 
*TDBAR,BBAR,PBAR 
1 FORMAT(10(/,21%,D15.9)) 


ore) 


convert time and intervals to radians 
TBARF=2.0*PI*FORBIT 
STEP=2.0*PI*STEP 
PRNT=2.0*PI*PRNT 


OO 


read user inputs of thrust angle and control parameters 
PRINT*, 'ENTER ALPHA’ 
READ*, ALPHA 
ALPHAR=ALPHA*PI/180.0D+0 


PRINT*, 'ENTER DLTRBR' 
READ*, DLTRBR 


PRINT*, 'ENTER ETHBR' 
READ*, ETHBR 


initialize counters 
TBAR=0.0D+00 
INDEX=0 
KOUNT=1 


AQ 


initialize state variables 
XBAR(1)=ROBAR 
XBAR(2)=0D+0 
XBAR(3)=0D+0 
XBAR(4)=(1/ROBAR)**1.5 
XBAR(5)=MOBAR 


ole) 


initially, set achieved orbit radii to center of orbital band 
RMNBR=ROBAR 
RMXBR=ROBAR 


set radius for thruster initiation 
RTHBR=ROBAR-DLTRBR 


CALCULATIONS 


QANAN AN 


cal] DRAG to determine initial and keplerian atmospheric drag 
CALLDRAG(BBAR,XBAR,BETABAR,VBAR,DBAR) 
DOBAR=DBAR 

C commence calculation loop. call DRAG to determine s/c velocity and atm drag 

100 CALLDRAG(BBAR,XBAR,BETABAR,VBAR,DBAR) 

C determine s/c energy 

EBAR=(VBAR*VBAR)Y/2.0-(1/XBAR(1)) 
C determine spacecraft thrust capability 
THBRCAP=TDBAR*DOBAR 
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OO O- -O! An os 


o) 


AQ 


AAF AY AM NAN AN 


determine keplerian thrust (=atmospheric drag) 
THBRKEP=DOBAR 
call CONTROL to fire or shutdown thrusters 
CALL CONTROL(ETHBR.XBAR,THBAR,RTHBR,EBAR.THBRCAP) 
call EQN to update the orbital equations of motion 
CALLEQN(DBAR,XBAR,XBDOT,VBAR,THBAR.ALPHAR,PBAR) 
call RK4 to do the runge-kutta integration over time 
CALL RK4(TBAR,XBAR,XBDOT,J,STEP,.INDEX) 


check for completion of runge-kutta integrator. if not done, loop back to 100 
IF (INDEX .NE. 0) GO TO 100 
runge-kutta complete: update variables for computation of orbital elements 
and output 
GAMMAR=ATAN(XBAR(2)((XBAR(1)*XBAR(4))) 
GAMMAD=GAMMAR*180.0D+0/PI 
compute flight path angle during thrusting, THGMA 
IF (THBAR .GT. 0.0D+00) THEN 
THGMA =GAMMAD 
ELSE 


THGMA=0.0D+00 
ENDIF 


RBAR=XBAR(1) 

MBAR=XBAR(5) 
MBRFUEL=MOBAR-MBAR 
FUELKEP=THBRKEP*PBAR*STEP 
TOTKEP=TOTKEP+FUELKEP 

IF (RBAR .LE. RMNBR) RMNBR=RBAR 
IF (RBAR .GE. RMXBR) RMXBR=RBAR 
ACTBAND=RMXBR-RMNBR 


call ORBELMTS to determine the classical orbital elements 
CALLORBELMTS(EBAR,RBAR, VBAR,GAMMAR,ECC,ABAR,APOGEE, 
*PERIGEE,PRDBAR) 


check if print interval acheived. if not, skip ouput section 
IF (KOUNT .LT. DNINT(PRNT/STEP)) GO TO 200 


OUTPUT 


WRITE(20,2)TBAR/(2.0*PI), RBAR,VBAR,MBRFUEL,EBAR 
FORMAT(2X,F6.2,2X,F11.9,2X,F11.9,2X,F11.9,2X,F1 1.9) 


WRITE(303)TBAR/(2.0*PI), ABAR,ECC,APOGEE, PERIGEE,PRDBAR 
FORMAT(2X,F6.2,1X,F11.9,3X,F8.6.3X,F11.9,1X,F11.9,2X,F8.4) 


WRITE(40,4)TBAR/(2.0* PI), DBAR,GAMMAD,THBAR,TOTKEP 
FORMAT(1X,F6.2,2X,F12.9,1X,F10.4,1X,F11.9,1X%,15.12) 


once printing completed, reset counter 
KOUNT=0 
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C update counter 

200 KOUNT=KOUNT+1 

Cc 

C check if simulation complete. if not loop back to 100 
IF (TBAR .LT. TBARF) GO TO 100 

Cc 


PRINT*, ‘Actual band achieved: ' 
PRINT*, ACTBAND 
Cc 


END 
CCCCCCCCCCCCCCCCCCCCCC 
C DRAG calculates the spacecraft velocity and the atmospheric drag acting on the 
C _ spacecraft, using a local exponential atmospheric density model 
Cc 


SUBROUTINE DRAG(BBAR,XBAR,BETABAR, VBAR,DBAR) 
IMPLICIT REAL*8(A-H,M-Z) 
DIMENSION XBAR(S5) 


calculate spacecraft velocity 
VBAR=((XBAR(2)*XBAR(2))+(XBAR(1)*XBAR(4))**2)**0.5 

calculate nondimensionalized atmospheric density 
RHOBAR=EXP(-BETABAR*(XBAR(1)-1)) 

calculate nondimensionalized atmospheric drag 
DBAR=RHOBAR*VBAR*VBAR/(2.0*BBAR) 


RETURN 

END 
CCCCCCCCCCCCCCCCCCCCCCO 
C EQN updates the spacecraft equations of motion 
Cc 


Yao OO 


SUBROUTINE EQN(DBAR,XBAR,XBDOT, VBAR,THBAR,ALPHAR,PBAR) 
IMPLICIT REAL*8(A-H,M-Z) 
DIMENSION XBAR(S5) .XBDOT(5) 


PI=[DATAN(1.0D+00)*4.0D+00 


A=XBAR(1)*XBAR(4)*XBAR(4) 
B=1/(XBAR(1)*XBAR(1)) 
C=DBAR*XBAR(2)/((XBAR(5)*VBAR) 
E=THBAR*SIN(ALPHAR)/XBAR(5) 


F=2.0*XBAR(4)*XBAR(2)/XBAR(1) 
G=DBAR*XBAR(4)/(XBAR(5)*VBAR) 
H=THBAR*COS(ALPHAR)(XBAR(5)*XBAR(1)) 


XBDOT(1)=XBAR(2) 
XBDOT(2)=A-B-C+E 
XBDOT(3)=XBAR(4) 
XBDOT(4)=-F-G+H 
XBDOT(5)=-PBAR*THBAR 
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RETURN 
END 
CCCCCCCCCCCCECCCCCCCCCC 
C ORBELMTS calculates the nondimensionalized orbital elements 
Cc 
SUBROUTINE ORBELMTS(EBAR,RBAR, VBAR,GAMMAR,ECC.ABAR, 
*APOGEE,PERIGEE,PRDBAR) 
IMPLICIT REAL*8(A-H,M-Z) 


Pl=DATAN(1.0D+00)*4.0D+00 


calculate angular momentum 
HBAR=RBAR*VBAR*COS(GAMMAR) 
calculateeccentricity 
PROB=(1.0+2.0*EBAR*HBAR*HBAR) 
IF (DABS(PROB) .LT. 1.0D-12) THEN 
ECC=0.0D+00 


Qa an 4 


E 
ECC=PROB**0.5 

ENDIF 
calculate semi-major axis 

ABAR=-1.0/(2.0*EBAR) 
calculate apogee 

POGEE=ABAR*(1.0+ECC) 

calculate perigee 

PERI erie 1.0-ECC) 
calculate period 

PRDBAR=2.0*PI*(ABAR** 1.5) 


RETURN 


END 
ececececececceccceccccce 
C —_RK4 does the fourth order runge-kutta integration 


Oe oO) Bk oO A 


Cc 
SUBROUTINE RK4(TBAR,XBAR,XBDOT,J,STEP,INDEX) 
IMPLICIT REAL*8(A-H,M-Z) 
INTEGER INDEX, ] 
DIMENSION XBAR(5),XBDOT(S5),SA VED(5),SAVEX(5) 
Cc 
INDEX=INDEX+1 
GOTO (1,23,4), INDEX 
Cc 


] DO 10 1=1J 
SAVEX(D=XBAR(I) 
SAVED(I)=XBDOT(I) 
10 XBAR(I)=SA VEX(1)+0.5D+0*STEP*XBDOT(I) 
TBAR=TBAR+0.5D+0*STEP 
RETURN 


DO 20 I=1J 
SAVED(ID=SA VED(1)+2.0D+0*X BDOT(I) 


NO 
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20 XBAR(D=SAVEX(1)+0.5D+0*STEP*XBDOT(I) 
RETURN 


Cc 
3 DO 30 I=1,J 
SAVED(D=SA VED(D+2.0D+0*XBDOT(D 
30 XBAR(D=SA VEX(1)+STEP*XBDOT(I) 
TBAR=TBAR+0.5D+00*STEP 
RETURN 


Cc 

4 DO 40 t=1J 

40 XBAR()=SAVEX(1D+STEP/6.0D+0*(SA VED(I)+XBDOT(I)) 
INDEX=0 
RETURN 


END 
CCCCCCCCCCCCCCCCCCCCCC 
C CONTROL determines if the spacecraft is within radial and energy limits prescribed 
C by the desired radial band and controls firing of the thrusters accordingly 
Cc 
SUBROUTINE CONTROL(ETHBR,XBAR,THBAR,RTHBR,EBAR, 
*THBRCAP) 
IMPLICIT REAL*8(A-H,M-Z) 
DIMENSION XBAR(5),XBDOT(S5) 


check if thrusters already on 
IF (THBAR .EQ. THBRCAP) GO TO 100 
check if orbit radius is less than or equal to thrusting radius 
IF (XBAR(i) .LE. RTHBR) THEN 
check if s/c energy is less than or equal to thrusting energy 
IF (EBAR .LE. ETHBR) THEN 
if both are true, then fire thrusters 
THBAR=THBRCAP 


Oreo CY 3 oeYQ 


C if not, do not fire thrusters 
THBAR=0D+0 
ENDIF 
E 


THBAR=0D+0 
ENDIF 
Cc 
C _ check if s/c energy is greater than thrusting energy. if so, turn off thrusters 
Me IF (EBAR .GE. ETHBR) THBAR=0D+0 


RETURN 
END 


Cc 
cccecceececccccccccecccccccecccceccececcccccecccccecccc 


Pea INPUT FILE 


D+00 
; -000000000D+00 
1.000000000D+02 
2.000000000D-04 
1.000000000D-01 
1.407286740D+02 
7.500000000D+01 
2.403901513D+04 
7.902309754D-01 
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ROBAR 
MOBAR 
FORBIT 

STEP 
PRNT 
BETABAR 
TDBAR 
BBAR 
PBAR 


B. SAMPLE OUTPUT FILES (through 20 orbits) 


SPACECRAFT CONFIGURATION: 234 
T/D:75 B:150 Isp: 1000 
REQUIRED BAND: 0.006035 
DLTRBR: 0.001035681 
ETHBR: -0.499482695 
ACHIEVED BAND: 0.0059946 


FILE: ORBPAR 
ORBIT__RADIUS 


0.999733807 
0.999457389 
0.999169910 
1.003238879 
1.003357471 
1.003534507 
1.003680096 
1.003794299 
1.003874691 
1.003919779 
1.003930433 
1.003904591 
1.003843996 
1.003747894 
1.003617000 
1.003452807 
1.003257664 
1.003034742 
1.002785952 
1.002513219 


VELOCITY __ FUEL BURNED 
0.000000000 


1.000133083 
1.00027 1285 


0.997016952 
0.997209856 
0.997430684 
0.997677081 
0.997947 130 


4l 


0.000000000 
0.000000000 
0.002423820 
0.002644308 
0.002913851 
0.003 182619 
0.003451645 
0.003720414 
0.003988924 
0.004257950 
0.004526460 
0.004795487 
0.005063997 
0.005332249 
0.005600501 
0.00586901 | 
0.006 138037 
0.006406547 
0.00667454 I 


ENERGY 


-.499522178 
-.499523400 
-.499524904 
-.499526248 
-.499527664 
-.499529162 
~.499530408 
-.49953 1863 
-.499533080 
-.499534523 
-.499536089 
~.499537673 
-.499539165 
~.499540455 
-.499543844 
-.499543844 


FILE: ORBEL 


ORBIT SEM-MAJ AX _ ECCTCTY _ APOGEE PERIGEE __ PERTOD 
1.00 0.999733728 0.000002 0.999736079 0.999731376 6.2807 
2.00 0.999457127 0.000005 0.999461975 0.999452278 6.2781 
3.00 0.999169364 0.000008 0.999176871 0.999161858 6.2754 
4.00 1.000956558 0.003012 1.004017094 0.997986175 6.2926 
5.00 1.000956558 0.003013  1.003972874 0.997940243 6.2922 
6.00 1.000954109 0.003008 1.003965144 0.997943075 6.2922 
7.00 1.000951095 0.003003 1.003957068 0.997945122 6.2922 
8.00 1.000948404 0.002999 1.003950002 0.997946805 6.2921 
9.00 1.000945566 0.002995 1.003943187 0.997947945 6.2921 
10.00 1.000942563 0.002991  1.003936143 0.997948982 6.2921 
11.00 1.000940067 0.002988  1.003930476 0.997949658 6.2920 
12.00 1.000937151 0.002985  1.003924458 0.997949844 6.2920 
13.00 1.000934713 0.002982  1.003919724 0.997949702 6.2920 
14.00 1.00093 1822 0.002980 1.003914742 0.997948902 6.2920 
15.00 1.000928684 0.002978  1.003909285 0.997948082 6.2919 
16.00 1.000925510 0.002975 1.003903572 0.997947449 6.2919 
17.00 1.000922520 0.002973 =1.003898246 0.997946794 6.2919 
18.00 1,000919936 0.002972  1.003894303 0.997945570 6.2919 
19.00 1.000916822 0.002971 1.003890236 0.997943408 6.2918 
20.00 1,000913145 0.002969 =: 1.003884904 0.997941385 6.2918 
FILE: ORBETC 
ORBIT___ DRAG GAMMA ___ THRUST FKT_ FUEL 
1.00 0.000021599 -0.0001 0.000000000 = 0.000103273109 
2.00 0.000022462 -0.0003 0.000000000 0.000206546218 
3.00 0.000023396  -0.0004 0.000000000  0.000309819327 
4.00 0.000013114 0.1157 0.000000000 0.000413092436 
5.00 0.000012893 0.1045 0.000000000 0.000516365544 
6.00 0.000012571 0.0888 0.000000000  0.000619638653 
7.00 0.000012313 0.0721 0.000000000 0.000722911762 
8.00 0.000012114 0.0546 0.000000000 0.000826184871 
9.00 0.000011975 0.0365 0.000000000  0.000929457980 
10.00 0.000011899 0.0179 0.000000000  0.001032731089 
11,00 0.000011880 -0.0009 0.000000000 0.001136004198 
12.00 0.000011924 -0.0197 0.000000000  0.001239277307 
13.00 0.000012028  -0.0382 0.000000000 0.001342550415 
14,00 0.000012194  -0.0563 0.000000000  0.001445823524 
15.00 0.000012424  -0.0737 0.000000000  0.001549096633 
16.00 0.000012718  -0.0902 0.000000000 0.001652369742 
17.00 0.000013078 -0.1056 0.000000000  0.001755642851 
18.00 0.000013500 -0.1197 0.000000000 —_0.001858915960 
19.00 0.000013988 -0.1324 0.000000000 0.001962189069 
20.00 0.000014543 -0.1434 0.000000000 0.002065462178 
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